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The spectrum of Gamma-ray Burst: a clue 
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ABSTRACT 

In this work we numerically calculate the thermal radiation efficiency of the baryonic outflow. 
The possible outflow acceleration in the transparent stage, which lowers thermal radiation 
efficiency, has been taken into account. In the standard internal shock model for the prompt 
emission, the fast shells should move with a typical Lorentz factor > 5T; otherwise the GRB 
efficiency will be in disagreement with the observations, where Tj is the bulk Lorentz factor 
of the shocked/emitting region. The photosphere radius of these fast shells is small and the 
thermal radiation is too strong to be effectively outshone by the internal shock emission. This 
is particularly the case for some extremely bright events having I 1 ; ~ 10 3 , like GRBs 0803 19B 
and 080916C. The absence of a distinct thermal component in the spectrum of most GRBs 
challenges the standard internal shock model and may suggest a non-baryonic (magnetic) 
outflow component. Though the magnetic outflow model seems favored by more and more 
data, it can hardly reproduce the typical GRB spectrum. In the photosphere-gradual magnetic 
dissipation scenario, the spectrum cuts off at ~ 1 GeV, too low to account for the observations 
of GRBs 0809 16C. In the sudden magnetic energy dissipation model, the low energy spectrum 
is expected to be F v oc j/ -1 / 2 , too soft to be consistent with the data F v oc We speculate 
that the low energy spectrum puzzle could be unveiled by the mechanism that particles, in the 
magnetic dissipation process, are repeatedly accelerated. 

Key words: hydrodynamics — gamma rays: bursts — radiation mechanism: nonthermal 



, ; 1 INTRODUCTION 



The mechanism that produces the prompt 7-ray emission in 
gamma-ray burst (GRBs) is still unclear. So is the physical com- 
position of the GRB outflows. In the standard internal shock sce- 
nario, the prompt soft 7— rays are the synchrotron radiation of 
the shock heated electrons and the outflows are baryonic. How- 
ever, there is an increasing interest in the magnetic fireball model, 
in which a considerable fract ion of the outflow energy is in the 
form of magnetic field (e.g.. |Usovlll992h iDuncan & Thompsonl 
1 1 9921 : lThompson1[T99 4 : Lyutikov & Bland ford 120031) . Quite a few 
pieces of independent evidences suggest that the GRB central en- 
gine might be strongly magnetized. First, the analysis of some 
well-studied optical flashes of GRBs reveal that the magnetic fields 
in the reverse-shock region are much stronger than that in the 
forward -shock region, so that the GRB outflows are probably mag- 
netized jFan et al. II2OO2I: Zhan g et al. I l2003l ;l Kumar & Panaitescu 
l2003tlGomboc et al. f2 009). Second, the absence of a distinct ther- 
mal spectrum component in most GRB s is consistent with the 
Poyn ti ng-flux dominated outflow m odel |]5aigne & Mochkovitch 
120021 ; iLvutikov & Blandford I l20olh . Third, the non-detection of 
bright optical flash in most GRB afterglows can be attributed to 
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a mild or high magnetiza t ion of the outflow dFan et al. I [20041 : 
IZhang & Kobavashi I l200l : iMimica et al. 1120091) . Fourth, the ab- 
sence of a GeV-TeV spectrum excess in the prompt emission of 
most GRBs detected by Fermi sate llite is in agreement with the 
magnetic fireball model iRmT2009). Last, the (possible) detection 
of the high linear polarization degree of some GRBs suggests that 
the magnetic fiel d involved in the synchrotron radiation could be 
globally ordered fcvutikov et aL~ll2003l : lGranotll2003l ; [ Gotz et al. I 
l2009h . 

Very recently, two discoveries rendered the magnetic fireball 
model more attractive. One is the detection and the successful op- 
tical polarimetry of the optical flash of GRB 090102. Its optical 
afterglow em ission declined as ~ i -1 6 and then got shallowed 
to ~ t~ 9 dGendre et al. 1 I2009T) . Such behaviors can be inter- 
preted as the weakly magnetized reverse-shock emission super- 
posed with the forward-shock emission. If correct, these optical 
flash photons would be moderately or even highly polarized. The 
ongoing polarizati on analysis seems to confirm such a speculation 
dSteele et al.ll2009l) . The data, however, has not been released yet, 
hampering us to go further. The other is the detection of the feature- 
less Band-type spectru m of GRB 0809 16 C in a very wide energy 
range 8 keV - 13 GeV dAbdo et alj2009h . If the GeV photons and 
the soft 7-rays were from the same region, the bulk Lorentz factor 
of the emitting region Ti has to be in order of 10 3 and the emit- 
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ting radius is R-, ~ 2r 2 c<5r. ~ 6 x 10 15 cm (Fi/10 3 ) 2 (<5r./0.1 s), 
where the typical variability timescale St is measured in the local 
frame of the burst dAbdo et al.ll2009l ; IZou et al. Il2009h and c is the 
speed of light. For suc h a large R-,, the widely discussed photo- 
sphere model of G RBs ^Thompson |l994l ; lRees & Meszaros 120051 ; 

er et al. I l2006h failed. The absence of a thermal component in 
the low energy band has been taken as a piece of evidence for the 
Poynting-flux dominated outflow model iZhang & Pe'er 2009), in 
which the initial radius of the outflow getting accelerated is taken 
as Ro ~ c8t ~ 10 9 cm. In principle, Ro can be as small as ~ 10 6 
cm, the physical size of a ste llar black hole or a magnetaf]. As re- 
vealed in lNakar et al. I d2005t) . the thermal radiation from a baryonic 
outflow depends on Ro sensitively. A small Ro can suppress the 
therm al emission effectively. With a similar argument, iToma et al. I 
(2009) suggested that a non-baryonic outflow component was not 
needed if Ro ~ 10 6 cm. In this work, we re-address that problem. 
We show that the thermal radiation efficiency of a b aryonic out- 
flow d oes i ncrease with Rp rapi dly, in agreement with Nak ar et al. I 
d2005h and IToma etail J2009I) . However, for GRB 080916C, as 
long as Ti ~ 10 3 , the standard internal shock model is hard to 
reproduce the data even for a Ro as small as ~ 10 6 cm. The phys- 
ical reason is the following (see section 2 for the detailed numer- 
ical approach). In the internal shock model, the fast shells carry- 
ing most of the energy should move with a bulk Lorentz factor 
Tf ~ 2r s hFi 3> Tj, where r s h is the strength of the internal shocks 
(plea se see sectio n 12.11 for the discussion). The photosphere ra- 
dius l|Paczvnskill990l : lDaigne & MochkovT tch 2002; Naka r et ail 
l2005h should be R ph ~ 6 x 10 20 cm (L/10 54 erg s _1 )r f - 3 , much 
smaller than R ph ~ 6 x 10 20 cm (L/10 54 erg S _1 )rr 3 that has 
been adopted in previous estimates, where L is the total luminos- 
ity of the baryonic outflow. It is well k nown that T{ ^ R p h/Ro 
jPiran et al" 1993; Meszaros et a l. 1 19931) , so we have 



T f < 5 X lOfLirR 



1/4 -1/4 



(1) 



Please note that the convention Q x = Q/10 x has been adopted 
in cgs units. Hence T sh » Tt/2T\ sC 2.5L 1 6 / 4 4 R ~^ /4 Tr£ t which 
is only marginally consistent with the GRB efficiency request (see 
footnote [2] for the details). 

This work is structured as the following. In section 2 we nu- 
merically calculate the thermal radiation efficiency of the baryonic 
outflow. We then discuss the implication of the observation of GRB 
080916C on the physical composition of its outflow. In section 3 
we discuss the spectrum problem in the magnetic fireball model 
and speculate about a possible solution. We summarize our results 
with some discussion in section 4. 



2 THERMAL RADIATION EXPECTED IN STANDARD 
INTERNAL SHOCK MODEL VS. THE DATA: 
SHEDDING LIGHT ON THE PHYSICAL 
COMPOSITION OF THE GRB OUTFLOW 

The thermal radiation from the GRB outflow has been widel 
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we focus on the thermal radiation efficiency. Our approach is as 
follows. We numerically solve the number, momentum and energy 
conservation laws of an extremely hot shell and get the evolution 
of the bulk Lorentz factor T, the comoving thermal energy density 
e', the comoving number density n and the observed surface tem- 
perature T bs = rr' (T' is the comoving surface temperature of 
the shell). The calculation stops when the radius reaches 

R ph w 6 x 10 18 cmiM^T" 2 , 

at which the thermal photons escape from the shell, where r\ = 
L/Mc 2 is the dimensionless entropy and M is the mass loading 
rate. For R > R p h (i.e., the transparent st age), the shell may sti ll 
be acceler ated by radiat i on v ia photon drag (Mes zaros et al. 1 1993b . 
Following iRossi et al. I u2006h . the acceleration of the outflow by 
the radiation in the transparent stage can be estimated as 



dT 



dR 167r7? 2 m p c 3 r 2 



O'T-fcph 



(1 



R 4 

R* ' 



(2) 



where the subscript ph represents the parameter at the photosphere 
radius, m p is the rest mass of the proton and <tt is the Thompson 
cross section. The final bulk Lorentz factor of the outflow is then 
given by 



1^3 ~ T^3 i 6Lph _2 
1 final ~ 1 ph H ~§T~ P 



(3) 



For Lth/L w e' ph /(e' ph + n' ph m p c 2 ) < 1 and then r ph ~ T], 
we have r nna i w F ph [l + 0Arie' ph /r ph (e' ph + n' ph m p c 2 )]. Hence 
the thermal radiation efficiency can be estimated by (usually the 
thermal radiation from the shell surface at R < R p h is ignorable) 



r?th « (1-0.4-H- 



In the case of L t h/L 



io^X 174 



r ph e' th + n' th m p c 2 

1 (i.e., i] > r ph 



(4) 



5 x 



), we have F nnal w r ph (l + 1.2^/rph) 1 / 3 



l.ir 2 n V /3 < r) and then t]th ~ 100%. 



1 If the rate of the accretion onto the nascent black hole is high up to ^ 
1 Mq s — 1 , the accretion disk flow becomes optically thick to neutrinos 
inside a radius ~ 10-R S , where _R S = 8.85 X 10 5 cm (Mbh/3Mq) is 
the black hole Schwarz child radius and Mbh is the mass of the black hole 
iPi Matteo et al. |2002|) . Most of the neutrino emission comes from outside 
this region. In such a case, we have Ro ~ 10-R S - On the other hand, the 
cooling of the disk material is dominated by neutrino radiation process, 
crucial for launching a baryonic outflow, only inside a radius ~ 10 8 cm 
iNaravan et al. Il200ll) . Therefore, in our calculation _Ro ranges from 10 6 
cm to 10 s cm, which is also consistent with what people find in the GR B 
spectrum modeling (e.g. jRvde et al. l200d|Pe'er |200alGao et al. |2009|) . 



2.1 The internal shock strength expected in typical GRBs 

For our purpose, the typical physical parameters of GRBs, in par- 
ticular the peak of the vF u spectrum e p , the 7— ray luminosity 
(L 7 ), the bulk Lorentz factor of the emitting region T; and the 
emitting radius 7? 7 , are needed. For the bright GRBs detected by 
BATS E, the distribution of e p peaks at ~ 200 keV jPreece et atl 
2000). For the bursts with a known redshift z detected so far, the 
averaged redshift is about 2. So in the burst frame, the typ ical 
peak e nergy should be (1 + z)e p ~ 600 keV. As shown in iLi I 
(2008), the distribution of L 7 for 64 long-duration Swift GRBs 
peaks at ~ 5 x 10 erg s~ . The bulk Lorentz factor has been de- 
rived in various ways (e.g..lLithwick & Sari Il200ll; l Molinari et al. I 
120071 ; IXue et al. II2OO9I ; IZou & Piranll2009l) and T; ~ 300 seems 
quite reasonable. _R 7 can be estimated by ~ 2Y 2 c5t. In this work 
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we adopt the intrinsic variability timescale of the prompt emis- 
sion St ~ 0. 05 s, as suggested by the num erical simulation of 
the collapsar dMacFadven & Wooslev I[l999h . So is expected 
to be in order of 10 14 cm. Larger is possible, as found in some 
previous estimates dZhang et al. l200q ; lLazz ati & Begelma n 1 20061 : 
iKumar et al.1l2007l) . 

Following Fan & Piran (2008), the magnetic field strength B 
and the typical random Lorentz factor 7 m can be estimated as the 
following. The 7— ray luminosity L 7 is related to the total luminos- 
ity of the emitting material L as L 7 ~ e e L/(l + Y), where e e (ee) 
is the fraction of the shock energy given to the electrons (magnetic 
field) and Y is the averaged Compton parameter. The comoving 
strength of the magnetic field can thus be estimated by 



B 



(2e B L/R;rfc 



2 p 2 sl/2 



4.5 x 10 4 Gauss (— + Y) * L* S2 R^uT^.s- (5) 



The synchrotron radiation frequency of electrons with a typi- 
cal Lorentz factor 7 m is ~ (1 + z)e p /h, which in turn yields 



7* 



and then 



r (l + z)e P 
[ 2.8 x W^hTiB 1 



11/2 



2000 



^(1 



&i(^)i[(i + y)L 7 , S2 ] 

600 kev £b 



'R. 



7,14' 



(6) 



[(p-l)/3(p-2)](e e /0.3)-\ 



7,14 

(7) 



where h is the Planck's constant, < £ c ^ 1 is the fraction of the 
electrons getting accelerated at the shock front, and p is the index 
of the power-law energy distribution of the accelerated electrons. 

One can see that even for £ c ~ 0.1, the shocks are relativis- 
tic, i.e., T s h ~ 2. Such an estimate is likely conservative since a 
larger f e ~ 1 is needed to acco unt for the early X-ray afterglow 
observations dFan & Piran Il2006l section 3.2 therein). Please note 
that r a h ^ 2 is also required to get a GRB efficienc)0 r/i ~ 20%. 
Therefore, in the internal shock model, for typical bright GRBs, we 
have Tf ~ 2r s hTi ~ 1000ri,2.5. For extremely bright bursts, like 
GRB 0803 19B and GRB 080 916C, r f ~ 5 x 10 3 is needed sinc e 
Ti ~ 10 3 and R 1 > 10 15 cm dZou et al. ll2009l : lAbdo et al.f2 009). 



so the total energy emitted during the acceleration phase is 



L thiS dR/(r 2 c) ~ 47vaR 2 T^(R /c) 



(9) 



where the relation T b s ~ const. dPiran et al. I 1 19931 : 
Mesza ros et al. 1 1 19931) has been taken into account, To is the 
temperature of the initial outflow and a is the Stefan-Boltzmann 
constant. 

The total energy of the shell can be estimated as 



Et t ~ AnRlcaTo St, 



(10) 



(ID 



where a = 4a /c is the radiation constant. 

The GRB efficiency contributed by the thermal emission leak- 
ing from the surface at R ^ R* is then given by 

£-th,s ^0 

E tot ~ AcSt 

Usually the central engine has a radius Ro ^ 10 6 cm. The typi- 
cal variability timescale of the GRB outflow may be mainly gov- 
erned by the accretion process and can be as long as ~ 50 ms 
dMacFadven & Wooslev |[l999h . For these typical parameters, we 
find 



TJth,. 



1.5 x 10~ 4 i? ,6(5t/50ms) _1 . 



(12) 



Therefore the thermal radiation of the accelerating shell is unim- 
portant unless St ~ Ro/coi R p h < R t . 

For fl« < R < R vh T' w r)~ Y T {R/ R f )~ 2,i 
(Pira n et al. I ll993l :l Meszar os et al. I ll993h . the subsequent thermal 
emission should have a luminosity Lth,s-l ~ 4tvi~i 2 R 2 aT' 4 ~ 
<±-kRIoT${R/R«)~ 2/3 and the detected emission should have a to- 
tal energy £ th , s -i = /£ ph 4 : TvR 2 aT^(R/R,)- 2/3 dR/(2 V 2 c) = 
3(R ph /R,) 1/3 R44-KRiaT^)/(2 V 2 c). Since (R ph /R,) 1/3 ~ 
10 -C T], we have iJth.s-i -C E t h,s- 



2.3 The thermal radiation efficiency of the baryonic outflow 

Following IFiran et al. I dl993h and lKonigl & Granotl d2002l) . for an 
extremely hot outflow we have the following number, energy and 
momentum conservation laws 



1 djn'T) 1 d(R 2 n'u) 
c dt + R? dR 



(13) 



2.2 The thermal radiation leaking from the surface at 

R < i? p h 

The acceleration of one baryonic shell with a width cSt (measured 
by the observer) is driven by the thermal photons an d can be ap- 
proximated as r ~ R/Rp fo r R < R* = T]Ro < i? p h dPiran et alTl 
Il993l : iMeszaros et al. |[l993h . The thermal emission from the sur- 
face of a shell in the case of R < R*(< R p h) can be estimated 
as 



47r.R 2 crT' 4 r 2 , 



(8) 



2 Let's just consider the collision of two shells. The masses and the Lorentz 
factors of the fast and slow shells are denoted as [Mi, M B ) and (I^f , r s ), 
respectively. The internal shocks are most efficient when an inner engine 
produces shells with compara ble energy but v ery different Lorentz factors, 
i.e., VfMf = r a M a iKobavashi et al. II 19971) . In such a case the merged 
shell has a bulk Lorentz factor T; ~ %/2r s and the efficiency is rg 
1 - (M f + M s )ry(r f M f + r s M s ) » 1 - {M ( /M s + 1)/V5 JPiranl 
1 19991) . Setting rj; ~ 20%, we have M { /M B ~ 0.14, T f ~ 7r a ~ 5T; and 

r sh - 2.5. 



1 d(wT 2 

c dt 



+ 



1 d[R 2 wTu] 



R 2 



OR 



icy 

c dt 



FA' 

c 



1 d(w'Tu) 1 d{R 2 w'u 2 ) _ dp' uM 



+ 



(14) 



(15) 



c dt R 2 dR dR c 

where the comoving entropy density w' = n'm p c 2 +7e', the ther- 
mal pressure denisty p' = (7 — l)e', u — Tf3 = \/T 2 — 1, and 
A' ~ aT' 4 /cTSt ~ e'/(4T8t). The specific heat ratio can be esti- 
mated by 7 w 1 + (aT' 4 /3 + n'kT')/(aT' 4 + 3n'kT'/2). In the 
current case we find that aT' 4 /3 3> n'kT', for which 7 w 4/3. 

In the case of F ^> 1, the above equations can be significantly 
simplified. Elimination of the r adiative cooling term from equations 
j 14b and J 1 5 b leads to (see also lKonigl & Granot 2002) 



dp' . r 2 a 1 df3 



where the convective derivative is d/dt = d/dt + f3cd/dR. 
Eqs. dl3t and d 1 5 1 > can also be approximated as 

1 d(Tn') 2Tn 



c dt 



+ 



R 







(16) 



(17) 
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middle and thin lines have 
been multiplied by a factor of 
(3, 1, 1/3), respectively. 
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Figure 1. The bulk Lorentz factor of the shell as a function of 
R. The solid, dashed and dash dot-dotted lines are for Ro = 
(10 6 , 10 7 , 10 8 ), respectively. The thick, middle and thin lines are for 
L = (10 54 , 10 53 , 10 52 ) erg/s, respectively, r) = 1000 is assumed in the 
calculation. 




The thick, middle and thin lines have 
been multiplied by a factor of 
(3, 1, 1/3), respectively. 



10 



10 



R (cm) 



Figure 3. The parameter e'/ (n'm p c 2 ) of the shell as a function of R.. The 
line styles and r) are the same as those of Fig[T] One can see that the larger 
the Ro , the higher the fraction of thermal energy escaping from the shell at 
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Figure 2. The observed surface temperature of the shell as a function of R. 
The line styles and T] are the same as those of Fig[T] 



and 



1 d{T 2 f3w') 
c 



dt 



+ 



2F 2 /3w' 
R 



utt 



(18) 



respectively. 

Eqs. d 161 i~8l l. together with the relation dR = /3cdt, are com- 
plete for solving T, e and n as functions of R. The starting point 
in our calculation is R — 5Ro, at which we take V — 5. The calcu- 
lation ends at R = R p h- 

Our numerical results have been plotted in Fig[T]-Fig[3] We 



find that T oc R, T ohs 



const, and e'/n'r 



oc R- 1 



for e' ^ n'm p c 



while T ~ const., T obs oc R~ 2/3 and 
e' /n'rripC 2 oc R ~ 2 ^ 3 for e' < n'm p c 2 . All are consistent with 
IPiran et ahl dl993l) . as expected. If at late times 7 approaches a con- 
stant lying between 4/3 and 5/3, with the general relation p' oc n' 1 
we have 



e jn oc R 



-2(7-1) 



oc R 



- 7 /2 



(19) 




log(v) 

Figure 4. A schematic plot of hiding the thermal emission by the non- 
thermal radiation powered by late energy dissipation, for example, internal 
shocks. 



2.4 Theoretical prediction versus the data: Constraint on the 
nature of the outflow 

Can the thermal emission be outshone by the nonthermal emission 
powered by internal shocks? For simplicity we assume the non- 
thermal emission component takes the form F v — Fo(y Ivq)~® 1 
for v < vo = e v /h and F v = Fa(v / Vo)~ for v > vq. 
As found in the data analysis, for typical GRBs f3\ ~ and 
f3 h ~ 1.25 jPreece et al. I feoOO). The thermal radiation peaks at 
a frequency ftfth, P = 2.82fcT Q b s and the corresponding flux is 
F„ th p = 2frz/ t 3 hjP /c 2 /[exp(frz/ t ii,p/fcTobs) - 1]. In the case of 
Pi < 1 and /3h > 1 (see Fig|4]for the details), in order to hide the 
thermal emission component, the ratio between the thermal emis- 
sion energy E t h and nonthermal emission energy E nt h should sat- 
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isfy 



£nth . 0.6(/3 h - A) 



E th " (1- A) (A - 1) 



^ "th.p ' 
1, 

/ ,y th 7 p ' 
V vn ' 



for v th>p < i/ ; 
for i^th.p = Vol 
for ^th, P > i*>. 

(20) 

For typical GRBs, E n th ^ 3 — 10£th is needed, otherwise the 
thermal component can not be hidden. The GRB internal shocks 
should have an efficiency r\i ~ 20% of converting the kinetic en- 
ergy of the outflow into radiat ion, as found in the afterglow mod- 
eling (e.g.. lFan&Piranll2006h . As a result, the thermal radiation 
efficiency should satisfy 77th ^ Tji/10 ~ 2%. Since the fast shells 
just take a fraction of the total energy^, the limit can be a bit higher. 
A reasonable requirement is 



'/th 



< 5%. 



(21) 



For L ~ 10 52 erg s~\ 77 ~ T f ~ 10 3 , and Ro ~ 10 6 cm, we have 
e ph /(e ph -|-n ph mpC 2 ) ~ 15% (see FigO, rj/T p h = 1.25, and then 
7/th ~ 7%, which violates the above request (though marginally). 
So the absence of a dis tinct thermal spectrum component in most 
GRBs (Rvde et al. ||20o3) may be a problem of the standard inter- 
nal shock model. The same con clusion has already been drawn by 
Daigne & Mochkovitch (2002). However in their modeling a very 
small £ c is needed and the resulting internal shock efficiency is very 
low a few percent). 

If GRB 080916C indeed had r f ~ 5 x 10 3 , the thermal ra- 
diation might be very strong. For this burst, the aftergl ow data had 
not b een collected until half a day after the trigger l lGreiner et al. I 
2009). The data are rare and can be well understood within the 
forward sh ock model supp osing the medium is a very weak stel- 
lar wind fcou et al. I 2009). The physical parameters can not be 
uniquely determined. The isotropic-equivalent kinetic energy of the 
outflow can be as large as ~ 4 x 10 55 erg and the corresponding 
GRB efficiency is r]; ~ 20% jGao et al. I2009T) . As shown in Figg] 
for L ~ 2 x 10 54 erg s" 1 and Ro ~ 10 6 cm, we have 77th ~ 10%, 
violating eq.J2U. A thermal spectrum component will be distinct. 
The absence of such a component m ay thus favor the non-b aryoinc 
(plausibly magnetic) outflow model. Izhang & Pe'erl ([2009) got the 
same conclusion. However in their approach Ro ~ 10 9 cm is as- 
sumed, much larger than what we adopt. 



3 THE SPECTRUM PROBLEM OF GRBS IN THE 
MAGNETIZED OUTFLOW MODEL 

The absence of a distinct thermal component is at odds with the 
standard internal shock model and may favor the magnetic fireball 
model. A self-consistent interpretation of the typical Band spec- 
trum of GRBs, however, is still unavailable, as shown below. 

In a pure photospher e model in the case of magnetar wind for 
GRBs, [Thompson I j 19941) showed that the typical Band spectrum 
of GRBs could be reproduced. However, for GRB 0809 16C, the 
prompt emission with a single power-law spectrum up to ~ 70(1 + 
z)~* GeV suggests that i? 7 ~ 10 16 cm, much larger than the site 
of th e photosphere ~ 10 9 cm suggested bv lThompsonl dl994l) . 

lGianniosl ( l2007l) showed that within the photosphere-gradual 
magnetic dissipation scenario, the low energy spectrum could be as 
hard as F„ oc u°. However, the corresponding high energy spec- 
trum was usually a bit harder than F v oc v~ x , inconsistent the data. 



3 For the most efficient internal shocks suggested in Kobayashi et al. 
1 19971) . V f M f ~ r s M s , so the fraction is ~ 1/2. 




Figure 5. The observed surface temperature of the shell (the thin lines) 
and e'/n'm p c 2 (the thick lines) as a function of R. The parameters are 
f] = 5000, Ro = 10 6 cm and L = 2 x 10 54 erg s" 1 . 



What's worse. lGiannios I < f2007t) predicted a cutoff at an energy ~ 1 
GeV, too low to account for the observations of GRB 0809 16C. 

Below we focus on the model of a sudde n magnetic energy 
dissip ation at Rj ~ 10 15 — 10 16 cm (e.g.. lLvutikov & Blandfordl 
2003). We introduce the parameter a to denote the ratio between 
the magnetic and the particle energy density. We perform a gen- 
eral study in which the details of the magnetic dissipation and the 
subsequent particle acceleration have been ignored. After the dissi- 
pation, the strength of the residual magnetic field can be estimated 
as 

B m ~ 300 Gauss fc^[a/(l + a)} 1 ' 2 L^R-^Y^ (22) 

where k 1 is the parameter reflecting the importance of the 
magnetic dissipation, which has been normalized to 0.1 because a 
residual magnetization o& ~ ka/[l + (1 — e e )(l — k)a] ^ 0.1 
may be needed to account for the absence of bright optical flashes 
in most GRB afterglows. Obviously k plays the same role of eb in 
estimating the strength of magnetic field of the emitting region. The 
electrons accelerated by the energy dissipation are assumed to take 
a power-law distribution oc y~ p for 7 e > 7 m - Similar to eq.© we 
have 



In 



2.5 x 10 4 [(1 + z)e p /600 keV] 



[*/(! + *)] 



^k~_l /4 L^ 



-l/4pl/2 



7,15.5- 



(23) 



The corresponding constraint on a reads 

a ~ 140Ce[3(p-2)/(p- l)]" 1 [(l + z)e p /600 keV] 1/2 

-V^fs/n , S M-1/4 D VS 

i.5- 



(e e /0.3)- 1 kZl /i L^[a/(l + a)]-^R!/ 2 



(24) 



The outflow has to be highly magnetized otherwise the dissipated 
energy is not enough to accelerate electrons to a typical random 
Lorentz factor ~ a fewx 10 4 . 

The cooling Lorentz factor can be estimated as 



7c 



50 L^[a/(l 



■a)] 1 fc_Ji? 7 ,i5.5rf i2 .B 



(25) 



We need a very large 
T; ~ 5000 Ll^la/il+otf^k^R 



-1/6 

'7,15.5 



(l+z)e p /600keV] 



1/6 



(26) 
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Figure 6. The synchrotron and synchrotron self-Compton spectrum of 
the electrons suffering significant inverse Compton cooling though within 
Klien-Nishina regime. The parameters are p = 2.5, j m = 7 X 10 3 , 
T; = 200, g ~ 35 and e c /k = 100. 



200 ll |2003l: IDerishev Il200l IWang et al. 1120091 ; iNakar et al. 1120091; 
Daigne et al. Il2009h . However, the resulting spectrum F v oc v 
for e c /k ~ 100 is not hard enough to account for the typical data 
F v oc u°. The lo w energy spectrum could be as hard as F v oc v~ 0A 
if e c /e B ~ 10 4 dNakar et al. Il2009b . corresponding to k ~ 10~ 5 . 
It is, however, unclear how the magnetic energy dissipation can be 
so efficient (i.e., k <C 1). It is also unclear whether the resulting 
spectrum for k ~ 10~ 5 can be well approximated by the Band 
function or not (see eq JA6t for the constraint). Observational tests 
of DKK's scenario may be available soon. 

A speculated solution to the low energy spectrum problem is 
the repeated acceleration of the particles in the energy dissipation 
process but much more work is needed to see whether it is the case. 
Some preliminar y discussion, but only f or unmagnetized outflow, 
can be found in IStern & Poutanenl fc004h . Please also bear in mind 
that in the multiple particle acceleration model the magnetization is 
required to be (much) higher than that suggested in eq.< !24t other- 
wise the re-accelerated electrons can not achieve a random Lorentz 
factor as large as ~ 10 4 (see eq.d23|l). 



to get 7 C ~ 7 m and then a low energy spectrum F v oc u 1 ^ 3 
that is roughly consistent with the observations. Such a large 
Ti is unrealistic (see eq.((TJ for the constraint) and is in con- 
tradiction with other constraints. So the low energy spectrum is 
likely F v oc v^ 1 ^ 2 , inconsistent with the data. Such an incon- 
sistence between the model and the data, already found in stan- 
dard internal shock m odel, is the so-called "low energy spectrum 
crisis of GRBs" (e.g., ICohen et al. 1 1 19971: iGhisellini et al. ||200C| ; 
iKumar & M cMahon 2008; |Piran et al. l2009h . ~ 

Eqs.( l24t and l |26t suggest a very low baryon loading of the 
outflow 

M jet ~ E jct /(rac 2 ) ~ 5 x 10" 10 M Ce" 1 r t ;2. B S je t,5 1) 

where Ej ct is the typical geometry-corrected energy of GRBs. It is 
unclear how such clean fireballs can be launched in the collapsar 
scenario. A hot massive neutron star as the GRB central engine is 
disfav ored because of the huge baryon pollution from such a star 
(e.g.. lLevinson & Eichler 1fl993h . 

In the above estimate of the synchrotron spectrum we have not 
taken into account the fact that the IC cooling of electrons is energy 
dependent. Such a correction may be crucial if the IC scattering 
process takes place in the Klein-Nishina regime. With Eq.l|23t we 
have g = 7 m (l + z)E- p /Tim c c 2 ^> 1. The factor g determines the 
regime of scattering of electrons with a random Lorentz factor 7 m 
on its synchrotron radiation photons at a typical energy e p . The IC 
cooling of the electrons with 7 C > rim c c 2 /(l + z)£ P is suppressed 
by the Klein-Nishina effect. As a result, the synchrotron emissiv- 
ity of these electrons increases with 7 e , which leads to a harder 
synchrotron spectrum and has been adopted to account for the typ- 
ical GRB X-ray spectrum F v oc u° (e.g.. IDerishev et al7ll200l[ 
hereafter DKK). However, as shown in the Appendix, the magnetic 
field in the emitting region should be very low otherwise the Klein- 
Nishina effect is too weak to modify the energy distribution of elec- 
trons and then the radiation spectrum (see eq.l |A6t ). Even for an 
unreasonably small k, say ~ 10~ 3 — 10~ 2 , DKK's scenario may 
still not work. In Fig|6]w e present our num erical spectrum based 
on the code developed in iFan et alH d2008l the instantaneous ap- 
proximation with some minor modifications). The l ow energy spec- 
trum does get hardened, as widely speculated (e.g. . IDerishev et al71 



4 CONCLUSION AND DISCUSSION 

In this work, we have calculated the thermal radiation efficiency 
of a baryonic shell (see section 2 for the details) and have dis- 
cussed the spectrum problem of GRBs (see section 3 for the de- 
tails). In the standard internal shock model for the prompt emis- 
sion, the fast shells should move with a typical Lorentz factor 
T f ~ 5ri ~ 10 3 ri, 2 .3 otherwise the GRB efficiency will be in dis- 
agreement with the observations. The photosphere radius of such 



fast shells is ~ 6 x 10 cm Z^r 



much smaller than that of 



the slow shells with a typical Lorentz factor < IV Consequently 
the thermal radiation from fast shells will be much stronger than 
that from the slow shells. In the internal shock model one should 
focus on the fast shells when investigating the thermal emission 
of GRB outflow. We find out that though most of the thermal en- 
ergy has been converted into the kinetic energy of the batyons, the 
residual thermal photons escaping from the surface at R ^ i? p h 
can not be ignored. The possible acceleration of the outflow in the 
transparent stage via photon drag may be able to lower the ther- 
mal radiation efficiency by a factor of ~ 0.5. For typical GRB pa- 
rameters (L, Ro, T f ) ~ (10 52 ergs" 1 10 6 cm, 10 3 ), we have 
a thermal radiation efficiency ~ 7% (see Tab[T]for a summary). 
These thermal photons are detectable and play an important role in 
cooling the electrons accelerated in the internal shocks. The non- 
detection of such a spectrum component in most GRBs thus chal- 
lenges the standard internal shock model. This is particularly the 
case for some extremely bright bursts with a featureless Band spec- 
trum and a very large Ti ~ 10 3 , like GRB 0809 16C. Please note 
that our conclusion is for Ro ~ 10 6 cm, the lowest value it could 
be. For a magnetized outflow, the ther mal emission is expected to 
be weak dLvutikov & Blandford Il2003h . consistent with the data. 

There is an increasing interest in the magnetic fireball model. 
A self-consistent interpretation of the typical Band spectrum of 
GRBs, however, is still unavailable (see section 3 for the details). 
In the photosphere-gradual magnetic dissipation scenario, the re- 
sulting spectrum cuts off at an energy ~ a few GeV, at odds with 
the observations of GRBs 080916C. In the sudden magnetic en- 
ergy dissipation model, the low energy spectrum is expected to be 
F v oc v~ x / 2 , too soft to be consistent with the data. It is sug- 
gested that the synchrotron radiation spectrum of electrons suffer- 
ing significant IC cooling but within Klein-Nishina regime can be 
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Table 1. The thermal radiation efficiency of fast shells involved in the standard internal shock model. 



typical Ti 


coiTesponding Tf 


Vth 


Vth 






for Ro = 10 6 cm 


for Ro = 10 7 cm 


bright GRBs (L 52 = 1) ~ 200 - 300 


~ 10 3 


~ 7% 


~ 18% 


extremely bright GRBs (L54 = 1) ~ 10 3 


~ 5 x 10 3 


~ 15% 


~ 35% 



much harder than the standard fast-cooling spectrum F v oc v ' 
jDerishev et al. | |2001 ). helping us solve the so-called low energy 
spectrum crisis of GRBs. However, to reproduce the data, the mag- 
netic field in the emitting region is required to be extremely low 
(see eq. JA6b and the last paragraph of section 3 for the discussion), 
which seems unrealistic in the magnetic dissipation scenario. The 
particle re-acceleration may be able to give rise to a harder low en- 
ergy spectrum but much more work is needed to see whether it is 
the case. 
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APPENDIX A: DKK'S SCENARIO AS A POSSIBLE 
SOLUTION OF THE LOW ENERGY SPECTRUM CRISIS 
OF GRBS: SOME REQUESTS 

It is widely believed that the prompt soft 7— rays are the syn- 
chrotron radiation of the shock heated electrons. With reasonable 
parameters, the synchrotron radiation peaks in soft gamma-ray 
band, consistent with the data. The resulting low energy X-ray 
spectrum F v oc v' 1 ^ 2 , is, however softer th an the typical low en- 
ergy spectrum F u oc u° JPreece et al. II2OO0T) . This inconsistency 
between the model and the da ta is sometimes called "low energy 
spect rum crisis of GRBs" (e.g-. ICohen et al. ll997l : lGhisellini et al. I 
l2OO0l) . 

In DKK's scenario (see Fig lAll for illustration), to obtain 



a hard X-ray spectrum F v 



rather than oc v x l 2 for 



£1 < hv < e p , the emitting electrons should satisfy: (1) 7c < 
7m(£p/£i)~ 1/ ' 2 < 7c < 7m- (2) The SSC cooling of the elec- 
trons at an energy ~ 7 m m c c 2 is important (the Compton parame- 
ter Y > 1) though in the Klein-Nishina regime. Let's show how it 
works. In general, th e energy distributi on of fast cooling electrons 
can be estimated by dNakar et al. I F2009) 



^oc[i + y( 7 c 



7. 



-(p+1) 



for 7c < 7c < 7m 
for 7c < 7m < 7c 



(Al) 



The synchrotron radiation of electrons with a random Lorentz fac- 
tor satisfying 7 C < 7c < 7m and Y(j e ) oc 7~ a > 1 has a spec- 
trum 



F v oc v 



(o-l)/2 



(A2) 



As long as the IC cooling of these "low" energy electrons is dom- 
inated by the photons in the spectrum segment F v oc i^ a ~ x ^' 2 , 
requiring that rim c c 2 /7 c ^ (1 + z)e p , we have 



y(7c) oc 



rim c 2 /h 7c F v dv 
Ub 



oc 7, 



-(o+l)/2 



(A3) 



where Ub is the magnetic energy density of the emitting region. 
Combing with the initial assu mption Y(y e ) oc 7e~ a > w e have a = 1 
and then reproduce eqs. (22) of lDerishev et al. I J200ll) 



dn/d r y c oc 7c 1 , F v oc v° . 



(A4) 



With eq. dAU . it is straightforward to show that in the energy 
range minfl/pK ( 7m)] 2 , [^(7m)] 2 }e P < hv < e p (see also 
iNakar et al. 2009) 

F v oc 



l/ -(P-i)/a j for y( 7m ) > i ; 
v- 1 ' 2 , for Y( 7m ) < 1. 



(A5) 



Therefore only for Y(jm) ~ 1 the X-ray spectrum can be 
approximated by F v oc u°. Y(7 m ) ^ 1 is also needed to satisfy 
the high energy radiation limit set by the recent Fermi observations 



L°9(yJ 



< rm c ly 




-Log(hv) 



Figure Al. The schematic plot of Derishev et al. (2001)'s model. 



(see lFan ]Hooi and the references therein). That's why below we 
focus on the ideal case Y(7 m ) ~ 1. In reality, the Band function 
jBand etal.lll993h is smooth across e p and 1/2 < Y(f m ) < 2 
is allowed by the data, which helps but does not solve all the 
fine-tuning problem. 

The request that Y(7 m ) ~ 1 imposes a tight constraint on the 
magnetization of the emitting region. For simplicity we consider 
that pi = and /3 h = p/2 ~ 1.25. The IC cooling of the elec- 
trons is dominated by the photons with an energy ^ rim c c 2 /7 c . 
The corresponding energy density of these seed photons is Uic « 



v°dv ~ FoTim e c 2 /h-y e . For7 e ~ j m , the 



request (c) reads that Uic ~ (ep/ei) 1 ^ 2 ^, with which we have 

F Ep/h w Jorim c c 2 //i7c,m ~ Ub(e v /e\) 1/2 , 

where the request rim c c 2 /7c,m ~ e P has been taken into account. 
The magnetic energy density Ub = B 2 /8tv and U 1 are related by 
Ub = £b(1 + Y)U~f/e c . The energy density of the soft gamma- 
ray emission can also be estimated by U~, ~ [p/(p — 2)]Foe p /h. 
Combining all these relations we have 



6B 



S)i(i + y). 



Usually E\ is below the low energy threshold of the detector ~ 10 
keV. The model thus demands 



±1 ~ 10o[_ J 

€b L 5(p-2) Jv 100ei 



]( _fp_)i ( i±r). 



(A6) 



This is a very general argument. If the low energy spectru m F v oc 
i/° is indeed produced in the way o ADerishev et al. I &200l\l . we can 
estimate e e /eb f equally e e /k) reliably without the need of any other 
information expect the prompt emission spectrum. 

With a reasonable e c ^ 0.3, eq,l lA6t gives that eb ~ 10~ 3 , 
which is reasonable in t h e ba ryonic outflow model (see 
IZhang. MacFadyen & Wang] 120091, for a recent simulation) 
but not for the magnetic fireball model. 

For completeness, below we present the other prediction of 
DKK's model, which can be tested observationally in the future. 
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Since we find out that DKK's scenario may not apply to the mag- 
netic fireball, here we turn back to the standard internal shock 
model. With eq.©, we have 



9 



r(l+*) 



p ( -)4 [( i + F)i 7 , 52 ] 



^7,14^1,2.5- 



600 keV 

Substituting eq.l lA6t into eq.© and eqs. JA7l l. we have 

(l + z)e Pl i r p 



(A7) 



7m ~ 5000 | 



600keV J L 5(p-2) J v 100ei 



-)°LJ 2 R 



"7,14' 

(A8) 



22 



. (l + z)e p a, 
600 keV 1 1 



2) J v 100ei 



l S ^7,52-R7,14 r i,2.5- 



(A9) 



One can see that for typical GRB parameters, the IC cooling of 
the electrons with energy ~ 7 m m c c 2 is indeed within Klein- 
Nishina regime (i.e., g 3> 1). On the other hand, the request that 
Ep/ei > g ^ (e p /ei) 1//2 should be satisfied otherwise the single 
power spectrum F v oc v° for ej < hv < e p can not hold. With 
eq.l |A9l >, such a request is equivalent to 

St 

J "f,52\ 



> 



130 



keV( 



-1/4 r 



-1/6 r 1/6 / 



-1/3 



1 + z """ v 100ei L 5(p-2) J 7 ' MV 0.05s' 

(A10) 

Therefore, if DKK's scenario is the solution of the low energy 
spectrum problem of bright GRBs, the X-ray flashes and X-ray 
flares should have a spectrum F v oc u~ 1 ' 2 rather than F v oc v°, 
which can be tested directly. The problem is that for such soft 
events , the low energy spect rum usually can not be reliably mea- 
sured dSakamoto et al. 112005b . 



